Low-disorder microwave cavity lattices for quantum simulation with photons 
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We assess experimentally the suitability of coupled transmission line resonators for studies of 
quantum phase transitions of light. We have measured devices with low photon hopping rates 
t/27r = 0.8 MHz to quantify disorder in individual cavity frequencies. The observed disorder is 
consistent with small imperfections in fabrication. We studied the dependence of the disorder 
on transmission line geometry and used our results to fabricate devices with disorder less than two 
parts in 10^. The normal mode spectrum of devices with a high photon hopping rate t/27r = 31 MHz 
shows little effect of disorder, rendering resonator arrays a good backbone for the study of condensed 
matter physics with photons. 
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Solving seemingly simple problems in quantum me- 
chanics is a formidable task for even the most sophis- 
ticated classical computers. For this reason Feynman 
proposed using controlled quantum systems to simulate 
and study other quantum systems [ ]. The development 
of such quantum simulators [ ] has since been an active 
area of research in a number of physical systems, includ- 
ing ultracold atoms in traps and optical lattices 
trapped ions [ ], and Josephson junction arrays [ ]. 

An idea of growing interest is to use photons as parti- 
cles in a quantum simulator for non-equilibrium systems 
[7-13]. According to this idea, a photon lattice is cre- 
ated with an array of cavity quantum electrodynamics 
(cQED) elements, each consisting of a photonic cavity 
coupled strongly to a two level system, or qubit. Sys- 
tems consisting of up to three coupled cavities have been 
realized for quantum information processing [14-16], and 
early proposals consider using larger arrays as a possible 
quantum computing architecture [^'^-^^], but a lattice- 
based quantum simulator of this type has yet to be real- 
ized. 

In these lattices, photons can hop between neighboring 
cavities and experience an effective photon-photon inter- 
action within each cavity, mediated by the qubit. The 
superconducting circuit architecture is an attractive can- 
didate for realizing such lattices due to the flexibility af- 
forded by lithographic fabrication and the relative ease of 
attaining strong coupling [ ]. Such cQED lattices have 
been predicted to exhibit a wide variety of phenomena, 
including a superfluid-Mott insulator transition [ - ], 
macroscopic quantum self trapping [ ], and fractional 
quantum Hall physics [ ]. 

In order to jump-start the implementation of these 
lattice-based simulators, we have fabricated and charac- 
terized 25 arrays of cavities, with each cavity designed to 
be identical. (Devices discussed here do not include qubits 
yet.) In this letter, we focus on understanding and re- 
ducing uncontrolled disorder in arrays of resonators in a 
kagome geometry, the most natural two-dimensional lat- 
tice for such transmission line resonators. We find that 
disorder in the individual resonator frequencies mainly 
originates from variations in the kinetic inductance due 



to small changes in the transverse dimensions of each 
resonator. We reduce disorder to less than two parts in 
10^ with a suitable choice in the geometric layout of our 
transmission line resonators. 

The system of coupled cavities is described by the 
Hamiltonian: 
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where a\ , a^ are the photon creation and annihilation op- 
erators corresponding to resonator i in the array. We 
denote the frequency of resonator i and its shift due to 
random disorder by uOr and 5i , respectively. The hopping 
rate between resonators i and j is given by 



U 



: 2ZQCij{uJr + h){^r + ^j), 



(2) 



where Cij denotes the coupling capacitance between the 
cavities and Zq the characteristic impedance of the trans- 
mission line [22] . The array we study consists of twelve 
resonators coupled capacitively at their endpoints in a 
two-dimensional kagome geometry (Fig. 1). Photon hop- 
ping was achieved by coupling triples of resonators in the 
interior of the array with symmetric three-way capaci- 
tors. This coupling scheme naturally results in a kagome 
lattice. 

Uncontrolled disorder is a serious impediment to us- 
ing such arrays for quantum simulation. Here, we work 
to understand the origins and minimize the level of dis- 
order to mitigate this problem. In an ideal array with 
uniform resonator frequencies and hopping rate t, the 
spacing between normal mode frequencies scales linearly 
with the photon hopping rate, with a frequency separa- 
tion between the highest and lowest modes of (3 + >/5)t; 
which arises from diagonalizing the 12x12 matrix result- 
ing from Eq. (1). Assuming that disorder in coupling 
capacitances is negligible, we find that disorder in res- 
onator frequencies leads to shifts of the normal mode fre- 
quencies through the first term in Eq. (1) by an amount 
~ 5i^ whereas the additional disorder in t 
in corrections ~ Sit/uOr- Since ujr ^ t 
proximate the photon hopping rate to be uniform with 



^j only results 
^j , we can thus ap- 




FIG. 1. a) (Color online) Device picture of twelve capaci- 
tively coupled resonators. The overlaid orange dashed lines 
have been drawn between the coupled resonators and illus- 
trate how the photonic lattice sites form a single kagome 
star. Transmission was measured between the ports labeled 
"Input" and " Output". b,c) Images of symmetric 3- way ca- 
pacitors with low hopping rate (t/27r = 0.8 MHz) with 10 /xm 
and 40 /xm wide center pins, d) Capacitor with high hopping 
rate (t/27r = 31 MHz) and 40 /xm wide center pin. e) Image of 
outer coupling capacitor [hi/2ti — 0.05 MHz) for 40 /xm center 
pin. f ) Cross-section of coplanar waveguide resonator with 
center pin width a, on a dielectric substrate Sr. 



value t = 2ZqC uj^ for nearest-neighbor resonators, and 
tij = for other resonator pairs. Therefore, we are pri- 
marily concerned with effects of disorder in resonator fre- 
quencies. 

Without disorder, there are eight distinct mode fre- 
quencies, four of which are doubly degenerate. The pres- 
ence of disorder breaks the degeneracy, widens the distri- 
bution of normal mode frequencies, and results in twelve 
distinct frequencies. We study the effects of disorder by 
numerically diagonalizing the Hamiltonian for random 
{Si} drawn from a Gaussian distribution with a stan- 
dard deviation a. The resulting histogram for the num- 
ber of eigenmodes N{u) du in a given frequency interval 
[cj, uj-\-duj] is shown in Fig. 2 for varying amounts of disor- 
der (J. When (J <^t^ disorder is negligible and the normal 
mode frequencies are all close to those of the ideal lattice. 
As a increases and becomes larger than t, the peaks in 
the distribution associated with individual normal mode 
frequencies broaden and ultimately merge. Once merging 
occurs, the observed mode frequencies and corresponding 
modes can no longer be easily identified with the idealized 
modes. In the limit of a ^ t, the normal mode histogram 
approaches a single Gaussian of width a from which the 
overall disorder of individual resonator frequencies can be 
extracted. For this reason, devices with a small hopping 
rate t are ideal for discerning the effects of disorder. 

We have fabricated and measured 25 arrays of twelve 
cavities to quantify disorder and assess the feasibility 
of quantum simulation in cQED lattices. By design, 
each coplanar waveguide resonator had a frequency of 
ujr/27r ^ 7 GHz, and an impedance Zq = 50 1^. The de- 
vices were fabricated using photolithography on 200 nm 
of Nb sputtered onto a 500 /im thick sapphire substrate. 
Each device was mounted using high performance silver 
paste and then wire bonded to a copper circuit board. 
Wire bonds were also used to connect all ground planes. 
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FIG. 2. (Color online) The normal mode histogram in the 
presence of disorder. Normal mode frequencies are calculated 
from Eq. (1) using a set of {5i} drawn from a Gaussian dis- 
tribution with standard deviation a. For each value of a, this 
procedure is repeated 10^ times. Histograms are generated 
from 10^ disorder realizations (for each value of a), and are 
normalized to the maximum number of counts for clarity. For 
a ^ t, the histogram is dominated by disorder and forms a 
single Gaussian. For a ^t, the histogram shows sharp peaks 
corresponding to the ideal normal mode frequencies. 



At the outer edges of the array, each cavity is capacitively 
coupled to a transmission line, resulting in a photon es- 
cape rate k = ^ZqC^^^uj^ to the continuum. This al- 
lows us to measure transmission through opposite ports 
(Fig. 1) of the array using a vector network analyzer. 
The unused ports were connected to 50 1^ terminators, 
though no significant difference was observed when the 
ports were left open. Each device was cooled to a base 
temperature of 20 mK inside a dilution refrigerator - a 
necessary requirement for future quantum simulations 
with small numbers of polaritons [ -11, 23]. 

The set of our 25 devices, summarized in Table I, in- 
cludes samples with two distinct hopping rates of t/27T = 
0.8 MHz and t/27r = 31 MHz. These nominal values were 
obtained from Eq. (2) by using values for the coupling 
capacitences determined using a finite-element calcula- 
tion. While the high-t devices allow us to access t ^ a 
and are most useful for quantum simulation, the low-t 
devices are the better choice for characterizing disorder. 

We extract normal mode frequencies from the peak 
positions in the measured transmission spectra (Fig. 3 
a)-c)) in order to determine the disorder. To account for 
small systematic shifts in devices made in separate fab- 
rication batches, all frequencies were expressed relative 
to the mean peak frequency of each spectrum. For low- 
t devices, not all twelve peaks are always visible. Such 
"missing" peaks can be due to normal mode degeneracies 
(occuring in the ideal case), as well as normal modes with 
small or vanishing amplitude in either of the resonators 
coupled to the input or output port. 



TABLE I. Results extracted from 25 measured devices. De- 
vices were characterized with two different photon hopping 
rates t and three different center pin widths a. The random 
disorder a was extracted from peak positions of the trans- 
mission spectrum for each device. The disorder decreases for 
increasing a. The ratio a/t is a metric of how the normal 
mode frequencies are effected by disorder. For the 40/im de- 
vices, a is reduced to less that two parts in 10^ of cJr/Syr. All 
uncertainties are computed from standard deviation of indi- 
vidual measurements. 

t/2n (MHz) a (/xm) cr/27r(MHz) a/t # Measured 
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FIG. 3. Transmission spectra for measured devices. The first 
column shows spectra for devices with a) t/27r = 0.8 MHz, 
a = 40 /xm, b) t/27r = 0.8 MHz, a = 20 /xm, c) and t/27r = 
0.8 MHz, a = 10 /im. The width of the spectrum decreases 
for increasing resonator width, demonstrating a decrease in 
a. The second column cl),e) shows transmission spectra for 
two nominally identical devices with t/27r = 31 MHz and 
a — 40 /im. Each scan contains twelve well defined peaks 
that are consistent between the two devices. Peak positions 
are similar to those expected, when accounting for a system- 
atic edge effect due to the difference between inner and outer 
capacitors. The inset shows the lowest energy mode that is 
localized on the inner six resonators in the absence of disorder. 



For low-t devices, analyzing the peak positions pro- 
vides a systematic method for extracting a from a trans- 
mission measurement. Specifically, the disorder strength 
can be extracted from the peak positions using: 




n\ 



dis\ 



1 ^ 

-^(n,-a)', (3) 



where n = 12 is the number of resonators in each sample, 
Qi and Qf^^ denote the twelve normal mode frequencies in 
the absence and presence of disorder, respectively. Cli and 



Qf^^ are their means (for a single disorder realization), 
whereas ensemble averages over disorder realization are 
denoted by (•). In the disorder less case, the "variance" 
of the normal mode frequencies of the kagome star is 3t^ . 

Applying this method to samples with a standard 
10 /im width of the transmission line center pin, we find 
that the disorder cr/27r = (9.1 ± 2.8) MHz is larger than 
expected from resonator length variations due to finite 
resolution in optical lithography. To investigate the ori- 
gin of this disorder, we fabricate devices with different 
widths a of the center pin, while maintaining a constant 
Zq throughout and find that there is a systematic depen- 
dence of disorder on a. 

The magnitude of disorder decreases with increasing 
center pin width (Fig. 4). This dependence of disorder 
on the device geometry can be attributed to random vari- 
ations in the width of the center pin that arise during mi- 
crofabrication. These variations in width result in varia- 
tions in the kinetic inductance L/^, which in turn affects 
the resonator frequency through the relation: 



1 



UJr = 



'^yij^m + Lk)Cu 



(4) 



where L^ is the intrinsic magnetic inductance and Ctot 
is the total capacitance. In normal metals, L^ is sup- 
pressed by electron scattering but in superconductors the 
DC electrical resistance is vanishing and L^ is no longer 
suppressed. Although L/^ is more relevant in supercon- 
ductors, it is still two orders of magnitude smaller than 
Lrm for the device geometry considered here. For a sin- 
gle resonator, Lk typically only results in a small shift in 
UJr [24, 25]. For arrays of coupled resonators, however, 
these small shifts can introduce significant disorder if the 
kinetic inductance contributions vary across the array. 

For the small length scales used here, the sensitivity of 
the kinetic inductance to variations in a decreases rapidly 
as the width a is increased [26] [27]. 

In our devices, we observe variations in the center pin 
width of up to ~ 600 nm and twice that for the dielectric 
gap, when examining them with a scanning electron mi- 
croscope. The random disorder expected due to kinetic 
inductance variations can be estimated by comparing uJr 
for cavities of equal length but with widths differing by 
the observed 600 nm, see Fig. 4. The random disorder 
observed here is consistent with variations in device ge- 
ometry and can be reduced to less than two parts in 10^ 
by making resonators with 40/im wide center pins. 

Using this strategy to reproducibly obtain devices with 
small disorder, we next turn to the high-t devices. Trans- 
mission spectra for all four of these devices revealed very 
similar normal mode frequencies, confirming that disor- 
der was small. Two representative transmission spectra 
are shown in Fig. 3 d),e). For all high-t devices, the 
lowest energy mode is significantly smaller in amplitude 
than the other eleven modes. In the absence of disor- 
der, the lowest energy mode is localized to the six inner 
resonators and cannot be driven from any port. For the 
infinite kagome lattice, it is these localized states that 
lead to the known flat bands [23, 28-30]. Disorder in the 




20 30 

Center pin width (|im) 



50 



FIG. 4. Random disorder versus center pin width for all 
devices. Disorder extracted from low-t devices is plotted in 
black with upward pointing triangles, while disorder extracted 
from high-t devices is plotted in green with a downward point- 
ing triangle. The curve shows the difference in frequency for 
two resonators, one with center pin width equal to the value 
on the horizontal axis and the other with a center pin width 
600 nm smaller and dielectric gap 1200 nm larger. Error bars 
are computed from standard deviation of individual measure- 
ments. 



array weakly breaks this localization and causes the mode 
to acquire a small amplitude in the outer resonators. 

For high-t devices, a is small compared to t and both 
variances on the right-hand side of Eq. (3) are large 
and nearly cancel each other. Consequently, an alter- 
nate method to extract a in these devices is used. In 
devices where t ^ cr^ the twelve peaks are easily identifi- 
able in the transmission spectra and directly indicate the 
variation in individual normal mode frequencies. In this 
limit, the frequency Q^^^ of the j^^ normal mode can be 
expanded to lowest order in the {Si} as 



n 



dis 



% 






(5) 



generate normal modes the two normal mode frequen- 
cies depend on mutually exclusive sets of the {Si}. Thus, 
these two eigenfrequencies fluctuate independently about 
a common mean value. In order to estimate disorder a 
from the measurements of the high-t devices, we first 
calculate the variance of the frequencies corresponding 
to each set of singly or doubly degenerate normal modes. 
Then, using Eq. (6), we scale this variance by the sum 
of the squares of the partial derivatives to calculate an 
estimate for the disorder a. Finally, we take the average 
of the estimates for a found for all of the sets of normal 
modes, weighted by the order of the degeneracy of each 
set. Using this method, we find a = (1.1 ± 0.6) MHz 
which is well into the limit of a <C t. 

The measured peak positions for the high-t devices 
agree well with expectations when accounting for a sys- 
tematic shift between inner and outer resonators due to 
the difference between inner and outer capacitances. In 
principle, the systematic shift can be compensated by 
adjusting the length of the edge resonators and will be 
of less importance in large arrays. As disorder is now 
understood, building larger arrays and adding qubits are 
the next logical and feasible steps. 

Photon lattices open the door for future experiments 
looking for quantum phase transitions and other many- 
body photon effects in coupled cQED arrays. Uncon- 
trolled disorder in arrays of bare cavities can be reduced 
to only a few parts in 10^. With random disorder suf- 
ficiently reduced, cavity arrays can be fabricated with 
controlled amounts of disorder in order to study local- 
ization effects. [31, 32]. While fabricating qubits on 
resonance with each cavity will be challenging, strong 
effective photon-photon interactions have already been 
observed even when qubits are off resonance [ ]. This 
strong off-resonant interaction should enable the possi- 
bility of studying many-body effects in the presence of 
disorder in qubit frequency, as resonant cQED is not a 
strict requirement. Therefore, quantum simulation with 
cQED lattices appears to be a realizable goal with super- 
conducting circuits. 



The variance of this normal mode frequency with respect 
to disorder is then 



ACKNOWLEDGMENTS 



Co 



{n[ 



E 



f) 






{nf) 



(6) 



The partial derivatives in Eq. (6) can be calculated nu- 
merically. Doing so, one finds that for each pair of de- 
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